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WASP- 10b: a 3Mj, gas-giant planet transiting a late- type 
K star 
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ABSTRACT 

>^ ■ We report the discovery of WASP-lOb, a new transiting extrasolar planet (ESP) 

discovered by the WASP Consortium and confirmed using NOT F1ES and SOPHIE 
radial velocity data. A 3.09 day period, 29 mmag transit depth, and 2.36 hour duration 
are derived for WASP- 10b using WASP and high precision photometric observations. 
Simultaneous fitting to the photometric and radial velocity data using a Markov-chain 
Monte Carlo procedure leads to a planet radius of 1.28Rj, a mass of 2.96Mj and ec- 
centricity of «0.06. WASP-lOb is one of the more massive transiting ESPs, and we 
compare its characteristics to the current sample of transiting ESP, where there is 
currently little information for masses greater than «2Mj and non-zero eccentrici- 
ties. WASP-10's host star, GSC 2752-00114 (USNO-B1.0 1214-0586164) is among the 
fainter stars in the WASP sample, with V=12.7 and a spectral type of K5. This result 
shows promise for future late-type dwarf star surveys. 

Key words: methods: data analysis - stars: planetary systems - techniques: radial 
velocities - techniques: photometric 
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Photometric transit observations of extrasolar planets (ESP) 
are important because the transit strongly constrains their 
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Figure 1. (a.) The top panel shows the SuperWASP light 
curve for WASP-lOb (1SWASP J231558.30+312746.4). All the 
data (apart from that from SuperWASP-N) were averaged in 
300 second bins. The data were phased using the ephemeris, 
T = 2454357.858081°' P = 3.09276 days. 



orbital inclination and allows accurate physical parameters 
for the planet to be derived. Their mass-radius relation al- 
lows us to probe their internal structure and is vital to our 
understanding of orbital migration and planetary formation. 
The radial velocity measurements which are used to confirm 
a candidate transiting ESP also provide more complete in- 
formation on the orbital eccentricity. 

As wide field photometric transit surveys have col- 
lected additional sky and tempor al coverage, and un- 
derstood their noi s e com ponents l)Collier Cameron et al.l 
l2007al : ISmith et all 120071 ) . the number of transiting ESP 
has g rown to ov er 50 in line with earlier predic- 
tions ijHornel 120031). Recent ly one such survey, Super- 
WASP ijPollacco et all l2006t ) published its first 5 con- 
firme d ESP, all of which have periods of less than 3 



days ( Andcrson "et alJ 



IPollacco et al 
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1.1 VZUU&). and reported 
120081 : iHebb et al.l |200S| ; 
IJoshi et al.ll2008l : IWest et al.ll2008n . SuperWASP is perform- 
ing a " shallow-but- wide" transit search, designed to find 
planets that are not only sufficiently bright (9 < V < 13) 
for high- precision radial velocity follow-up to be feasible on 
telescopes of moderate aperture, but also for detailed stud- 
ies such as transmission spectroscopy during transits. De- 
tails of the WA SP project and obse rvatory infrastructure 
are described in lPollacco et all (|2006h . 

In this paper we present the WASP photometry of 
1SWASP J231558.30+312746.4 (GSC 2752-00114), higher 
precision photometric follow-up observations with the MER- 
CATOR and Tenagra telescopes, and high precision radial 
velocity observations with the Nordic Optical Telescope new 
Fibre-fed Echelle Spectrograph (FIES) and the OHP SO- 
PHIE collaboration. These observations lead to the discov- 
ery and confirmation of a new, relatively high mass, gas- 
giant exoplanet, WASP-lOb. 



2 OBSERVATIONS 

1SWASP J231558.30+312746.4 (GSC 2752-00114) was mon- 
itored by SuperWASP-N starting in 2004. SuperWASP is a 
multi-camera telescope system with SuperWASP-North lo- 
cated in La Palma and consisting of 8 Canon 200-mm f/1.8 
lenses each coupled to e2v 2048 x 2048 pixel back illumi- 
nated CCDs. This combination of lens and camera yields a 



field of view of 7.8° x 7.8° with an angular size of 14.2" per 
pixel. During 2004, SuperWASP was run with 4 or 5 cam- 
eras as the operations moved from commissioning to routine 
automated observing. We show the WASP-lOb SuperWASP 
light curve in Figure 1. 



2.1 Higher Precision Photometry 

We obtained photometry of WASP- 10 with the MEROPE 
instrument on the 1.2m MERCATOR Telescope in V-band 
on 1 September 2007. Only a partial transit was observed 
due to uncertainties in the period and epoch from the Su- 
perWASP data. Observations were in the V-band with 2x2 
binning over ~2.9 hours. Despite clear conditions, exposure 
times were varied from 25-30s to account for changes in see- 
ing and to keep below the saturation limits of the chip. This 
allowed 170 images to be taken. There was a drift of only 
1 binned pixel in x and y on the chip during the run. The 
MEROPE images were first de-biased and flat fielded with 
combined twilight flats using IRAF and the apphot package 
to obtain aperture photometry of the target and 5 nearby 
companion stars using a 10 pixel radius. Finally, the light 
curve was extracted and normalized to reveal a depth of ~33 
mmag. 

Further observations of WASP-10 were taken as part of 
an observing program sponsored by the Las Cumbres Obser- 
vatory Global Telescope NetworlQ on the Tenagra II, 0.81m 
F7 Ritchey-Chretien telescope sited in the Sonora desert in 
S. Arizona, USA. The science camera contains a lk x lk 
SITe CCD with a pixel scale of 0.87 arcsec/pixel and a field 
of view of 14.8' x 14.8'. The filter set is the standard John- 
son/Cousins/Bessel UBVRI set and the data presented here 
have been taken in I band. 

Calibration frames were obtained automatically every 
twilight, and the data were de-biased and flat-fielded using 
the calibration section of the SuperWASP pipeline. Object 
detection and aperture pho tometry was then performed us- 
ing daophot (|Stetsonll2007l ') within IRAF. Differential pho- 
tometry was derived from a selection of typically 5 to 10 
comparison stars within the frame. 

These confirmed the object had a sharp egress with 
an amplitude of 0.033±0.001mag. The MERCATOR V and 
Tenagra I light curves show consistent transit depths, con- 
firming that the companion is a transiting ESP. 



2.2 Spectroscopic Follow-up 

We obtained high precision radial velocity (RV) follow-up 
observations of WASP-10 with the 2.5m Nordic Optical Tele- 
scope (NOT) new Fibre- fed Echelle Spectrograph (FIES), 
supplemented with observations from the Observatoire de 
Haute-Pr ovence's 1.93 m teles cope and the SOPHIE spec- 
trograph (jBouchv et al] [2006). We present a summary of 
the FIES and SOPHIE RV data in Tabled] 



2.2.1 NOT and FIES 

Spectroscopic observations were obtained using the new 
FIES spectrograph mounted on the NOT Telescope. A to- 
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Table 1. Journal of radial- velocity measurements for WASP- 
10 (1SWASP J231558.30+312746.4, USNO-B1.0 1214-0586164). 
Stellar coordinates are for the photometric apertures; the USNO- 
B1.0 number denotes the star for which the radial-velocity mea- 
surements were secured. The quoted uncertainties in the radial ve- 
locity errors include components due to photon noise (Scction l2.2jl 
and 10 ms _1 of jitter (Section 13.21 1 added in quadrature. 



D, J D 


^cxp 

( b ) 




Vi 
km s~ 


i 




FIES 














2454437.540 


2400 


-11. 


.028 


± 


0. 


.026 


2454463.377 


2400 


-11. 


.941 


± 





.030 


2454465.342 


2400 


-11. 


.003 


± 





.018 


2454466.335 


2400 


-11. 


.804 


± 





.021 


2454490.329 


2400 


-11. 


.013 


± 


0. 


.143 


2454490.358 


2400 


-10 


.990 


± 





.120 


2454491.340 


2400 


-11. 


.955 


± 





.024 


SOPHIE 














2454340.569 


3300 


-11. 


.657 


± 





.008 


2454342.505 


3600 


-11. 


.575 


± 





Oil 


2454508.262 


1680 


-11. 


.027 


± 





.014 


2454509.268 


1680 


-11. 


.336 


± 





.017 


2454510.276 


1680 


-11. 


.990 


± 


0. 


.020 


2454511.262 


1680 


-11. 


.135 


± 





.016 


2454512.262 


1680 


-11. 


.244 


± 





.016 



tal of seven radial velocity points were obtained during 2 
December 2007, 28-31 December 2007 and 24-25 January 
2008. WASP-10 required observations with an exposure time 
of 2400s due its relative faintness (V=12.7) yielding a peak 
signal-to-noise ratio per resolution element of «60-70 in the 
Ha region. FIES was used in medium resolution mode with 
R=46000 with simultaneous ThAr calibration. We used the 
bespoke data reduction package FIEStooJf] to extract the 
spectra and a specially developed IDL line-fitting code to 
obtain radial velocities with a precision of 15-25 m s _1 (ex- 
cept for the poor night of 24 January 2008, JD 2454490). 

2.2.2 OHP 1.9 m and SOPHIE 

Additional radial velocity measurements were taken for 
WASP-10 on 2007 August 29 and 30, and again between 
2008 Feb 11 and 15 wit h the OHP 1.93 m t elescope and the 
SOPHIE spectrograph l|Bouchv et al.ll2006t ). a total of 7 us- 
able spectra were acquired. We used SOPHIE in its high 
efficiency mode, acquiring simultaneous star and sky spec- 
tra through separate fibers with a resolution of R=40000. 
Thorium- Argon calibration images were taken at the start 
and end of each night, and at 2- to 3-hourly intervals 
throughout the night. The radial-velocity drift never ex- 
ceeded 2-3 ms _1 , even on a night-to-night basis. Although 
errors for each radial velocity measurement are limited by 
the photon-noise. Thus, the average radial velocity error is 
~ 14 ms -1 and includes the 2-3 ms _1 systematic error and 
the contribution from the photon-noise. Typical signal-to- 
noise ratio estimates for each spectra were «30 (near 5500 
A). The SOPHIE WASP-10 spectra were cross-correlated 
against a K5V template provided by the SOPHIE control 

3 http:/ /www. not. iac.es/instruments/fics/ficstool/FIEStool. html 



Table 2. Stellar parameters for WASP-10. The last 5 parameters 
were derived from the SME analysis of the FIES spectroscopy. 

Parameter WASP-10 



RA (J2000) 


23 15 58.3 


Dec (J2000) 


+31 27 46.4 


V 


12.7 


distance 


90±20 pc 




4675 ± 100 K 


logs 


4.40 ± 0.20 


[M/H] 


0.03 ±0.20 


v sini 


<6 kms -1 


"rad 


-11.44 ± 0.03 kms 



and reduction software. Typical FWHM and contrasts for 
these spectra were «10.2-10.4 km s _1 and 30-31%, respec- 
tively. The cross-correlation techniques and derivation of er- 
rors in the ra d ial ve locity measurements are presented in 
IPollacco et all l|2008h . 



3 RESULTS AND ANALYSIS 
3.1 Stellar parameters 

We merged all available WASP-10 FIES spectra into one 
high-quality spectrum in order to perform a detailed spectro- 
scopic analysis of the stellar atmospheric properties. Radial 
velocity signatures were carefully removed during the pro- 
cess. This merged spectrum was then continuum-normalized 
with a low order polynomial to retain the shape of the broad- 
est spectral features. The total signal-to-noise ratio of the 
combined spectrum was ~180 per pixel. We were not able to 
include the SOPHIE spectra in this analysis, because these 
spectra were obtained in the HE (high-efficiency) mode, 
which is known to suffer from problems with removal of the 
blaze function. 

As previously un dertaken for our analysis of WASP-1 
l|Stempels et al.ll2007l ). and WAS P-3 dPollacco et al.ll2008ri 
we employed the methodology of IValenti fc Fischerl ||200*5l ). 
using the same tools, techniques and model atmosphere 
grid. We used the package Spectroscopy Made Easy (sme) 
l|Valenti fc Piskunovl 1 19961 ) . which combines spectral syn- 
thesis with multidimensional \ 2 minimization to determine 
which atmospheric parameters best reproduce the observed 
spectrum of WASP-10 (effective temperature T c g, surface 
gravity logg, metallicity [M/H], projected radial velocity 
wsini, systemic radial velocity « ra d, microturbulence v m ic 
and the macroturbulence t) mac ). 

The four spectral regions we used in our analysis are: (1) 
5160-5190 A, covering the gravity-sensitive Mg b triplet (2) 
5850-5950 A, with the temperature and gravity-sensitive Na 
l D doublet; (3) 6000-6210 A, containing a wealth of differ- 
ent metal lines, providing leverage on the metallicity, and (4) 
6520-6600 A, covering the strongly temperature-sensitive H- 
alpha line. In addition we analyzed a small region around 
the Li I 6708 line to possibly derive a lithium abundance, but 
no Li I 6708 was detected for WASP-10. The parameters we 
obtained from this analysis are listed in Table [5] In addi- 
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tion to the spectral analysis, we also use available photom- 
etry (from NOMAD, TASS4 and CMC14 catalogues), plus 
2MASS to estimate t he effective temperature using the In- 
frared Flux Method (|Blackwell fc Shallidll977h . This yields 
T c ff = 4650 ± 120 K, which is in agreement with the spectro- 
scopic analysis and a spectral type of K5. The characteristics 
of WASP- 10 are also given in Tabled 



3.2 Markov-chain Monte Carlo analysis 

Transit timing and the radial- velocity measurements provide 
detailed information about the orbit. We modelled WASP- 
lOb's transit photometry and the reflex motion of the host 
star simultaneously using the Mar kov-chain Monte-Carlo 
algorith m described in detail by ICollier Cameron et al.l 
i|2007al ), an d the same technique s that were applied to 
WASP-3 by IPollacco et all (|2008h to which we refer the 
reader for more details. 

We find WASP-lOb to have a radius 1.28±g"_g|J2j, mass 
of 2.96±{J;^Mj and a significant non-zero eccentricity of 
0.059lg;°j4 . The best fit solution for the MCMC model for 
a circular orbit (e = 0) has a y 2 55 higher than the solution 
with non-zero eccentricity, and thus, the eccentricity is sig- 
nificant at > 99.6% confidence level using the F test. The 
values of the parameters of the optimal solution are given, 
together with their associated 1-a confidence intervals, in 
Table H The FIES+SOPHIE radial- velocity data measure- 
ments are plotted in Figure [5] together with the best-fitting 
global fit to the SuperWASP-N, MERCATOR, and Tenagra 
transit photometry. 



3.3 Line-bisector variation 

Line bisectors have been shown to be a powerful di- 
agnostic in distinguishing true extra-solar planets from 
blen ded and eclipsing s t ellar systems c hrom ospheric activ- 
ity l|Queloz et all l200ll ). iTorres et ail (|2004l ) showed, that 
for OGLE-TR-33 line asymmetries which changed with a 
1.95 day period, it was a blended system. From the cross- 
correlation function (CCF) we obtained the line bisectors 
and these are plotted, as a function of RV, in Figure [3] 

We quantified the significance of the bisector variation 
as follows. We determined the inverse-variance weighted av- 
erages of the RV and bisector span as 



b = 



where the and bi are the RV and span bisector values 
respectively and the weights Wi are the inverse variances of 
the individual bisector measurements. The uncertainty in 
the span bisector is assumed to be 2.5 times the uncertainty 
on the RV in our data. If we define Xi — Vi — v and yi — bi — b, 
then the slope is determined as 



Var(a) 



Ei 



The value of the scaling factor a is determined with signal- 
to-noise ratio 
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Figure 2. Simultaneous Markov-chain Monte Carlo (MCMC) 
solutions to the WASP-10 photometry and radial velocity data, 
a. Top panels show the MCMC solutions to the combined 
SuperWASP-N, MERCATOR V, and Tenagra / band photome- 
try, b. The lower panel shows the MCMC solution to the FIES + 
SOPHIE radial velocity (RV) data. (FIES RVs are shown as filled 
circles and SOPHIE as filled squares). The model fit to the RV 
data includes orbital eccentricity (solid line), and for a circular 
orbit (dashed line). Both RV models also include the Rossitcr- 
McLaughlin effect, which is small for this system given the low 
vsini of the host star (<6kms — 1 ). 
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Figure 3. Line bisectors as a function of radial velocity (RV) 
for WASP-10. Plot symbols are the same as Figure 2. Analysis 
of these line-bisectors for WASP-10 does not show a correlation 
between the bisector velocity (V spa n) and stellar radial- velocity 
(see text). 



We obtain SNR = 1.16, indicating a non-significant corre- 
lation between the bisector span and radial velocity varia- 
tions. This demonstrates that the cross-correlation function 
remains symmetric, and that the radial- velocity variations 
are not likely to be caused by line-of-site binarity or stellar 
activity and indicate WASP- 10b is an exoplanet. 



4 DISCUSSION 

Photometric surveys have now provided a large sample of 
transiting ESP that can be used to determine their mass- 
radius relation and provide constraints on their composi- 
tions. Here we presented the discovery of a new ESP with a 
mass of 2.96Mj, 1.28Rj radius, and a significant eccentricity 
of 0.059t°;ow- We now discuss the properties of WASP-lOb 
in relation to the current sample of transiting ESP, starting 
with its non-zero eccentricity. 

Most of the current sample of published transit- 
ing ESP have orbits consistent with being circular and 
are fit with models using zero eccentricity as is ex- 
pected for short-period planets in orbits with semi- 
major axes < 0.2 AU. Recent work jjackson et al.l l200Sl : 
iMardling 2007 and references therein! ) has investigated the 
effects of tidal dissipation on the orbits of short period ESP. 
The evolution of the orbital eccentricity appears to be driven 
primarily by tidal dissipation within the planet, giving a cir- 
cularisation timescale substantially less than 1 Gyr for typ- 
ical tidal dissipation parameter, Q p = 10 5 to 10 6 . WASP-10 
is a K dwarf with a spin period of 12 days and J— K=0.62 
and is rotating more slowly than star s of comparable colour 
in the Hyades dTerndrup et alj|2000h . This suggests a rota- 
tional age between 600 Myr and 1 Gyr. Thus, the persistence 
of substantial orbital eccentricity in WASP-lOb is therefore 
surprising. 



Table 3. WASP-10 system parameters and 1-<t error limits de- 
rived from MCMC analysis. 



Parameter 



Symbol 



Value 



Transit epoch (BJD) 
Orbital period 
Planet/star area ratio 
Transit duration 
Impact parameter 

Stellar reflex velocity 
Centre-of-mass velocity 
Orbital semimajor axis 

Orbital inclination 
Orb. eccentricity 
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Stellar mass 
Stellar radius 
Stellar surface gravity 
Stellar density 

Planet radius 
Planet mass 

Planetary surface gravity 
Planet density 
Planet temp (A = 0) 
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One plausible mechanism for maintaining the high ec- 
centricity is secular interaction with a n additional planet in 
the system. lAdams fc Laughlinl (|2006t ) explore the effects of 
dynamical interactions among planets in extrasolar plane- 
tary systems and conclude outer planets can cause the in- 
ner planet to move through a range of eccentricities over 
timescales that are short when compared to the lifetime of 
the system, but very long when compar ed to the current 
observ ational baseline. However, recently iMatsumura et all 
(2008) have argued that an unseen companion driving short- 
Period systems is unlikely. They present an upper limit of 
1 Mjveptune for a possible unseen companion in the GJ 436 
system and exclude this base d on the current radial v eloc- 
ity upper limits of ^ 5 m/s. Mat sumura et alJ (|2008h also 
present a range of tidal quality Q p timescales that could 
be as large as 10 9 years, and argue that this new class of 
eccentric, short period ESP are simply still in the process 
of circularizing. WASP-lOb has not been extensively stud- 
ied to rule out a putative outer plane that may be driving 
its eccentricity. Thus, the f»6% eccentricity of WASP-lOb 
makes it an attractive target for future transit-timing varia- 
tion studies, and for longer-term RV monitoring to establish 
the mass and period of the putative outer planet. 

The majority of transiting ESP found have masses be- 
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low 1.5Mj, although there are a few more massive ESP. HD 
17156, and COROT-Exo-2 have similar masses to WASP- 
10b and although there are two more massive ESP, the 
nearly 9 Mj HAT P-2 (HP 149026b1 ([Bakos et all 120071 ) 
and 7.3 Mj WASP- 14b (|joshi et al.ll200Sr ), this higher mass 
region has been poorly explored. Additional transiting ob- 
jects in the mass range are important for completing the 
current ESP mass-radius relations and constraining their 
compositions. The current sample of transiting extrasolar 
giant planets (ESP) reveals a large range of densities. We 
derive a mean density for WASP-lOb of wl.89 g cm' 3 (1.42 
pj) and it would lie along the higher density contour in a 
mass-radius plot l|Pollacco et al.ll2008l ; ISozzetti et al.ll2007t) . 

One ultimate goal of our transit-search programme is 
to provide the observational grist that will stimulate and 
advance refined models for the fo rmation and evolution of 
the hot and very hot Jupiters He.s. Burrows et all 1 19971 ; 
iFortnev et al ] 120071 : ISeaeer et alJl2007h . By thus constrain- 
ing the underlying physics, we will have a richer context for 
the interpretation of the lower mass planets expected from 
missions such as COROT and Kepler. 
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